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'■ Abstract 

In this review we show that the multi-particle scattering amplitudes in J\f — 4 SYM at large N c and 
in the multi-Regge kinematics for some physical regions have the high energy behavior appearing 
from the contribution of the Mandelstam cuts in the complex angular momentum plane of the 
corresponding t-channel partial waves. These Mandelstam cuts or Regge cuts are resulting from 
gluon composite states in the adjoint representation of the gauge group SU(N C ). In the leading 
logarithmic approximation (LLA) their contribution to the six point amplitude is in full agreement 
' with the known two-loop result. The Hamiltonian for the Mandelstam states constructed from n 

^ , gluons in LLA coincides with the local Hamiltonian of an integrable open spin chain. We construct 

• the corresponding wave functions using the integrals of motion and the Baxter-Sklyanin approach. 

This article is an invited review for a special issue of Journal of Physics A devoted to Scattering 
Amplitudes in Gauge Theories. 
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1 Introduction 
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At high energies s 3> — t in QCD the elastic scattering amplitude for the process AB — > A'B' in the 
leading logarithmic approximation (LLA) 



a s In s ~ 1 , a s <C 1 (1) 

has the Regge form [TJ 

A 2 ^2 = 2g8 XA x A ,T c AA — gT% B , 5 Xb x b , , t = -q 2 . (2) 

Here T c are the generators of the gauge group SU (N c ), X r are the particle helicities and j(t) = l+co(t) 
is the gluon Regge trajectory for the space-time dimension D = 4 — 2e 

In the framework of the dimensional regularization the parameter \x is the renormalization point for 
the 't Hooft coupling constant 

where 7 = — V'(l) is the Euler constant and ip{x) = (lnr(x))'. The gluon trajectory j(t) was calculated 
also in the next-to-leading approximation in QCD [2] and in the SUSY gauge models [3]. 
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In LLA gluons with momenta k r (r=l,..,n) are produced in the multi-Regge kinematics 

s = (pa+ Pb) 2 > s r = (kr + fc r -i) 2 > -U = <ir , k r = q r+1 - q r , (5) 
where the amplitude has the factorized form (see also section [2]) 

A 2 ^ 2+n = 2s8 XAXA ,gT2 Al ^^gC,(q 2 , qi )e;(k 1 )T c % i ^^...^±l gT%£ 5 XXb , . (6) 

1l <?2 1n+l 

Here C^(q 2 ,qi) is the effective Reggeon-Reggeon-gluon vertex. In the case when the polarization 
vector e fl (ki) describes the gluon with a positive helicity in its cm. system with the particle A' one 
can obtain [3] 

C=C,(q 2 ,q 1 )e;(k 1 ) = V2^, (7) 

where the complex notation q = q x + iqy for the two-dimensional transverse vector q was used. 

The elastic scattering amplitude with vacuum quantum numbers in the t-channel is calculated in 
terms of the production amplitude A 2 ^ 2+n with the use of the s-channel unitarity [TJ. In this approach 
the Pomeron appears as a composite state of two Reggeized gluons. It is convenient to present the gluon 
transverse coordinates in the complex form together with their canonically conjugated momenta [HE] 

* ■ . d * . d . . 

Pk = x k + iy k , Pk = x k - iy k , Pk = i-^— , Pk = ■ ( 8 ) 

CPk op k 

In this case the homogeneous Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation for the Pomeron wave 
function can be written as follows [1] 

E^(p u p 2 ) = H l2 y(p 1: p 2 ) , A = -^-2 min E, (9) 

where A is the Pomeron intercept entering in the asymptotic expression for the total cross-section 
at ~ s A . The BFKL Hamiltonian has a rather simple operator representation [5] 

#12 = In \p lP2 \ 2 + -^r(ln \pi 2 \ 2 )pip 2 + -^(ln |pi 2 | 2 )PiP2 - 4^(1) (10) 

with pi2 = P\ — P2- The kinetic energy is proportional to the sum of two gluon Regge trajectories 
u(~\pi\ 2 ) (i = 1,2). The potential energy ~ In \p\ 2 \ 2 is obtained by the Fourier transformation from 
the product of two gluon production vertices C^. This Hamiltonian is invariant under the Mdbius 
transformation [6] 

apk + b 

Pk -> —7 , 11 

cpk + a 

where a, b, c and d are complex parameters. The eigenvalues of the corresponding Casimir operators 
are expressed in terms of the conformal weights 

1 it 1 Tl 

m = — h iv H — , m = — V iv , (12) 

2 2 2 2 v ' 

where v and n are respectively real and integer numbers for the principal series of unitary representa- 
tions of the Mdbius group SL(2, C). The eigenvalues of H\ 2 depend on these parameters and can be 
written in the holomorphically separable form [8] 

E m ,m = £m + £m, £» = 1p(m) + Tp(l - 771 ) - 2^(1) , (13) 

where ip(x) = (lnT^))' . 
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The Pomeron intercept in LLA is positive 

Q> 

A = 4— iV c ln2>0 (14) 

7T 

and therefore the Froissart bound at < c In 2 s for the total cross-section is violated [1] . To restore the 
broken s-channel unitarity one should take into account the contributions of diagrams corresponding to 
the i-channel exchange of an arbitrary number of reggeized gluons in the i-channel. The wave function 
of the colorless state constructed from n reggeized gluons can be obtained in LLA as a solution of the 
Bartels-Kwiecinski-Praszalowicz (BKP) equation [7\ 

E^ = ¥ , A = min E. (15) 

Air 

In the N c — > oo limit the color structure is simplified and the corresponding Hamiltonian has the 
property of the holomorphic separability [8j 

n 

= ^ H k>k+1 = hf® + h<-°> , [h® , hP>] = . (16) 
k=l 

It is a consequence of the similar property for the pair BFKL hamiltonian Hyi (jlOj) and the energy 
E m ,fh CGI- 

The holomorphic Hamiltonian in the multi-color QCD can be written as follows (cf. (|10B ) 

^ (0) = E4°l+i . h< $ = HP1P2) + - (mpi 2 )pi + - (ln Pl2 )p2 ~ 2^(1) , (17) 
k Pl P2 

where ip(x) = (lnT(x))'. As a result, the wave function ^ has the holomorphic factorization [8j 

^> = J2 <h,r ^(Pl, -, Pn) ^(Pl, P* n ) , (18) 
r,r 

which in the case of two-dimensional conformal field theories is a consequence of the infinite dimen- 
sional Virasoro group. Moreover, the holomorphic hamiltonian is invariant under the duality 
transformation [9] 

Pi -)• p i<i+1 ->■ Pi+i , (19) 

combined with its transposition. 

Further, there are integrals of motion q r commuting among themselves and with UOj : 

4 0) = Pk 1 k 2 Pk 2 k z ---Pk 2 kf-Pk r k 1 Pk 1 Pk 2 ---Pk r , [q r ,h]=0 . (20) 

fcl<fc2<...<fc r 

The integrability of the BFKL dynamics in LLA and its relation with the Baxter XXX-model was 
established in ref. [10j . This remarkable property is related also to the fact that h coincides with 
the local Hamiltonian of an integrable Heisenberg spin model |11| (see also |12|). Eigenvalues and 
eigenfunctions of this hamiltonian were constructed in refs. |13|. [14] in the framework of the Baxter- 
Sklyanin approach [15] . 

In the next-to-leading approximation the integral kernel for the BFKL equation was calculated in 
Refs. [31 [16]. In QCD the eigenvalue of the kernel contains the Kronecker symbols <5 nj o and 5 n ,2 but 
in N = 4 SYM it is an analytic function of the conformal spin and has the property of the maximal 
transcendent ality [2l[T7]. This extended N = 4 super symmetric theory appears in the framework of 
the AdS/CFT correspondence [HI [191 [20] . It is important, that the eigenvalues of one-loop anomalous 
dimension matrix for twist- 2 operators in M = 4 SYM have the maximal transcendentality property 
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and are proportional to the expression —ip{j — 1), which is related to the integrability of evolution 
equations for quasi-partonic operators in this model 2 [j. The integrability persists also for some op- 
erators in QCD |22j. The maximal transcendentality principle suggested in ref. [17] gave a possibility 
to extract the universal anomalous dimension up to three loops in M = 4 SYM [231 [24"] from the cor- 
responding QCD results [25] . The integrability of the M = 4 model was demonstrated for anomalous 
dimensions of other operators in higher loops and at large coupling constants |26 [ I27 [ [28] . In partic- 
ular, the asymptotic Bethe ansatz allowed to calculate the anomalous dimensions in four loops [29J. 
This result is in an agreement with the next-to-leading BFKL predictions after taking into account 
the wrapping effects [30] . The maximal transcendentality was helpful for finding a closed integral 
equation for the cusp anomalous dimension in this model {311 [32] with the use of the 4- loop result [33] . 
The thermodinamic Bethe ansatz and the approach based on the Y-systems allows to calculate the 
spectrum of anomalous dimensions in M = 4 SYM for an arbitrary coupling constant [34] . 

In recent years a new line of investigations has been started, which also shows remarkable properties 
of N = 4 SYM: the study of scattering amplitudes. A few years ago, Bern, Dixon and Smirnov (BDS) 
suggested a simple ansatz for the gluon scattering amplitudes in this model [35]. This ansatz was 
verified for the elastic amplitude in the strong coupling regime using the AdS/CFT correspondence [36] . 
But the BDS hypothesis does not agree in this regime with the calculation of the multi-particle 
amplitude [37] . leading to the conclusion that a non- vanishing remainder function, R> n \ has to exist 
which provides the necessary corrections to the BDS amplitudes. The property of the conformal 
invariance of the BDS amplitudes in the momentum space was discussed in ref. [38], and the relation 
with the Wilson loop approach was suggested in ref. [39] generalizing the results of the strong coupling 
calculations of ref. [36]. However, in ref. [40] it was found that the BDS amplitudes A n for n > 6 in 
the multi-Regge kinematics do not have correct analytic properties compatible with the Steinmann 
relations [41| . It is a consequence of the fact, that these amplitudes do not contain the Mandelstam 
cuts [40] . The cut contribution was obtained from the BFKL-like equation for the amplitude with 
the f-channel exchange of two reggeized gluons in the adjoint representation of the gauge group [30] . 
This equation was solved in LLA and the two-loop expression for the 6-point scattering amplitude in 
the multi-Regge kinematics was derived [42] . The two- loop correction to the remainder function was 
calculated numerically for some values of external momenta in an agreement with expectations based 
on the Wilson loop approach [13]. In a recent paper [44] . by solving the set of Y-equations, also the 
strong coupling limit of the remainder function has been studied. 

The existence of the Mandelstam cut contribution, found first for the 6-point amplitude, gener- 
alizes to multiparticle amplitudes with n > 6. As an example, the 2 — > 2n amplitude will contain 
Mandelstam-cut contributions composed of n reggeized t-channel gluons. These n gluon t-channel 
states can be expressed in terms of solutions of the BKP-like equation in the adjoint representation. 
Most remarkable, in LLA the corresponding Hamiltonian is integrable: it coincides with the local 
Hamiltonian of an integrable open Heisenberg spin chain [45] . 

In this review we present a summary of the Mandelstam cut contributions to the inelastic scattering 
amplitudes in N = 4 SYM and their properties of integrability. We first review the analytic structure 
of n-point amplitudes in the multi-Regge kinematics and describe the main features of the Mandelstam 
cut contributions. In the subsequent section we compare our results in LLA with the exact two- loop 
calculations of Goncharov, Spradlin, Vergu and Volovich (GSVV) [46] and consider a relation with 
collinear kinematics. The rest of the review is devoted to the integrability of the BKP Hamiltonian 
in the adjoint representation. 

2 The analytic structure of scattering amplitudes in the Regge limit 

Let us begin with a brief summary of the analytic properties of scattering amplitudes in the multi- 
Regge limit. It is well-known that the 2 — > 3 amplitude in the multi-Regge kinematics with the 
exchanged reggeons having definite signatures Tj = ±1 in the crossing channels t\ and t<i can be 
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written as a sum of two terms 



pole 

siHisreufc < 2 2 c} 2 + n^m^iCi < 2 j 4 2 , «i2 = k 2 = ^ , (2i) 



1VI 2^3 _ |„ |wia|„|W2£ f ,W2„12 , L |W21U|Wit t ^,12 ^ _ 1,2 S 1 S 2 



r(ti)r(i 2 ) 

where 

^ = e -i«wi _ n ? ^ = g-fcrwa _ T2 > = e -i^ 12 + TiT2 ^ 6i = g-i^a + ^ _ (22) 

r(tj) are the residue functions of the exchanged Regge poles, and k a is the transverse momentum of 
the produced particle. In (|21|) we assumed that s,s±,S2 and k±2 are measured in some characteristic 
mass /i 2 . 

The gluon Regge trajectories in N = 4 SYM can be written as (see ([3D) 

/ 2x 7ir i 0? _ , 9 2 ^ c /OQ x 

Ui = w(-q- ) = — — In — , 7k 4a , a = — , w i2 = u>i - u> 2 , (23) 

4 A z 07T Z 

where ij = — q 2 , is the cusp anomalous dimension and A 2 ~ /i 2 exp(l/e) for D = 4 — 2e with 
e — > —0. The parameter A 2 can be considered as an effective mass of gluon. The real coefficients 
c\ 2 , c^ 2 in M = 4 SYM are obtained from the BDS amplitude and given below [ID] 

12 I xi I sinvr(a;i - u a ) 12 . , sinvr(a; 2 - u a ) , . 

Ci = 1 12 : , c 2 = li2 : , (24) 

Sin 7TW 12 Sin 7TW 21 

where the Reggeon-Reggeon-gluon vertex Ti 2 in the physical region s, si, s 2 > is 

k 2 A 2 

Ti 2 (ln Ku - in) = \T 12 \ exp(i7r u> a ) , ui a = — In -f-* , (25) 

ln]r 12 | = — (-^n ^—ln + - In — 1 n + ^ J • (26) 

The amplitude representation in (|21|) is compatible with the Steinmann relations |41j forbidding 
the simultaneous singularities in the overlapping channels s± and s 2 . The expression in ()21[) can be 
rewritten in the factorized form 

M 2^l _ Liuit. T/T1T2 |„ _|W 2 £_ T/T1T2 _ ^12 12 , &1 



.... , .sir&^hrfc, ^ = ^c| 2 + ^ c ^. (27) 

1 (Cl)l (12) £l ?2 

The BDS-amplitude which holds in the planar approximation can also be written as a sum of two 
terms: 

M BDS 

r(f \rl S = (-«i)" 12 (-^i2) W2 cl 2 + (- S2 ) aJ21 (-^i2) Wl 4 2 , (28) 

where we put the normalization point fj, 2 in the Regge factors equal to unity. 

For the 2—7-4 amplitude in the multi-Regge kinematics the situation is more complicated. In 
agreement with the Steinmann relations, the Regge-pole scattering amplitude can, again, be written 
as a sum of 5 terms as illustrated in Fig. [TJ 

For the signatured amplitude, this representation is equivalent to the factorized form: 

= \ si rciv nv2 | S2 r 2 6^ r2T3 is 3 re 3 , (29) 

where V T2T3 is obtained from V T1T2 (|27p with the corresponding substitutions 

F r 2 r 3 = |3 c 23 + |2 c 23 (3Q) 
?2 £3 
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Figure 1: The analytic representation of the 2—7-4 scattering amplitude. The dashed lines denote 
possible energy discontinuities. 



For the second produced gluon with the transverse momentum k\, the coefficients c 23 and phase 
read 



r 23 



23| 



sin7r(o;2 — u>b) 

Sin 7TW23 



r 23 
c 2 



23| 



sm-7r(a;3 — ojf, ) 
sin 7ra;32 



7A' , k 2 A 2 2 

° q2^3 



S2«3 



S123 



(31) 



In the planar approximation one expects that the scattering amplitude, M^^, in accordance with the 
Steinmann relations [5T], has the form [40] : 



M 



pole 



2^4 



(_ S1 )«" (- S0 12^12r 3 (-^12^2 3 )" 3 Ci 2 cf 



r(ti)r(t 3 ) 

+ (_ S3 )-3 2( _ Sl23K23) - 21 (- SKl2K23 )^ c f c f + (- S K 12 « 23 )^ ("Sin 2 ("S3)" 32 Ci 2 C 23 

C 2 C x 



, / wit \wz t smTTWi sm7rw 2 3 12 23 

+ {S2) (-S012K12) 13 (-SK12K23) " " 



+ (-S 2 )" 23 (-S123K2 3 )" 31 (-SK 12 «23)' 



sin7TW2 sm7ro;i3 

sin7ro;3 shim^i 12 23 

— : c 2 c i 

sm7ro;2 sm7ra;3i 



(32) 



A closer look at the last two terms shows that a model with Regge-poles only exhibits unphysical 
poles, indicating that a pure Regge model maybe incompatible with the correct analytic structure 
of multiparticle amplitudes in the multi- Regge kinematics. In fact, the LLA analysis of N=4 SYM 
gauge theory has shown that, in addition to the gluon Regge pole, there exists also a Mandelstam 
cut in the complex angular momentum plane, which removes this inconsistency. The Mandelstam 



Figure 2: Contributions of a) the Regge poles and b) the Mandelstam cut to the 2 — > 4 scattering 
amplitude in the t2-channel. The wavy lines represent reggeized gluons. The Mandelstam cut appears 
bound state of two reggeized gluons. 

cut appears in the angular momentum plane of the ti channel and describes a bound state of two 
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or more reggeized gluons as depicted in Fig. [3 In the planar approximation, it shows up in the 
special physical kinematic regions where the invariants in the direct channels have the following signs 
s, s 2 > 0; si, s 3 , s i2> $123 < or s,s 1 ,s 2 ,s 3 < 0; s 012 , si 23 > [ID], and it is not visible in the physical 
kinematic region where all energies are positive. In the following these special physical regions will be 
named "Mandelstam regions". 

The 2—7-4 amplitude in the multi-Regge kinematics can be written as a sum of the Regge pole 
and Mandelstam cut contributions [40] 

M 2 ^ 4 = M|!| e 4 + M 2 c ^ 4 , (33) 

where is a generalization of two last terms in (|32p , and it is non-zero only in the two kinematic 

regions restricted by the inequalities s,S2 > 0; si, S3, S012, S123 < and s, s%,s 2 ,S3 < 0; S0127S123 > 0. 
There is some freedom in redistributing terms between the Regge pole and the Mandelstam 

cut M|^ 4 contributions. Using this fact and representation (132f) one can write in the region s, s 2 > 
0; Sl,S 3 ,S 12,Sl23 < [17] 



M. 



pole 

2— »4 —iivu)2 



|siM*2M*3Mri2||r23|r(ti)r(i 3 ) 

and 



e~ 1 ™ 2 cosiru ab (34) 



' - ' ' ~ i e~ inuJ2 / — f(u 2 >) e~ in ^' \s 2 \ u * , (35) 



|ai| wl |s2| < *|a3| W8 |ri2||r23|r(ti)r(t 3 ) J- loo 2m 

where co a b is obtained from (p5|) and ([3T]) and reads 



IK k|q| 
k 2 n 2 



u a b = — In^-f-f. (36) 



In the other physical region, where s, si,S2,s 3 < 0; soi2,si23 > (corresponding to the physical 
channel for the 3 — > 3 transition) we find |47] 

ole 

cos iruj a b (37) 



M pole 



Iai| wi |s2| wa |a3| ws |ri2||r23|r(t 1 )r(t 3 ) 

and 

M|^ 4 . r°° duj 



Lw*'^' (38) 



|8i| w i|s2| wa |s3| W8 |ri2||r23|r(ti)r(t 3 ) 

The function f{oo 2 >) is pure real and denotes the partial wave in the complex angular momentum 
plane. 

The origin of this 'restricted' appearance of the Mandelstam cut contribution M|!^ 4 can be traced 
back to Mandelstam's argument for the existence of the Regge cuts (Fig. [3]): If we put for the reggeon 
momentum k = apA + @Pb + k±, it is easy to see that, in the planar (large- _/V c ) limit with all energies 
being positive, the integrations over a and (3 have singularities only in the upper half planes and lead 
to the absence of the Mandelstam cut contribution. However, if by pulling the produced particles to 
the left we 'twist' the reggeons (ladders) in the t\ and t% channels (Fig. [4]), there will be singularities 
on both sides of the real a and f3 axis, and the Mandelstam singularity remains. Note that, despite 
this 'twisting', the amplitude is still planar. Returning to the sum of the five contributions in Fig. [H 
it can be shown that the Mandelstam cut contribution should be present only in the last two terms: 
taking into account the phase structure one finds that, in the physical kinematic region, where all 
energies are positive, the cut cancels in the sum of the two terms. This cancellation does not work, if 
we are in the 'mixed' regions s, s 2 > 0; s±, s 3 , S012, S123 < or s, si, s 2 , s 3 < 0; S012, S123 > 0. 
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Figure 3: The diagrammatic structure of the Mandelstam cut. 




Figure 4: Twisting a planar 2—^4 amplitude. 



Turning now to the BDS amplitudes, this Mandelstam cut contribution is missing and must there- 
fore be contained in the remainder function, i?2-i>4- It is believed that the full MHV amplitude in the 
planar (large- N c ) approximation can be written in the factorized form: 

M 2 _ >4 = Mf4f i? 2 ^ 4 . (39) 

Indeed, in [42] it was shown that in the region s,S2 > 0; s%, S3, S012, S123 < the correct form of the 
LLA 2—7-4 scattering amplitude is 

M^t = M 2 B 4f (1 + iA^i), (40) 

where M^f is the BDS amplitude [35] and 

AS _ . £ ( _ ir r _^ my-* c^y +5 l^,, _ a . (41) 

Here k a , /c^ are the complex transverse components of the produced gluon momenta, q±,q2, <?3 are the 
momenta of reggeons in the corresponding crossing channels, and 

(Tl Tl \ 

2-0(1)- -0(1 + iv + ~) - + iv - ~)) (42) 

is the eigenvalue of the BFKL Hamiltonian in the adjoint representation. The correction A^i^ is 
Mobius invariant in the momentum space, and it is important to note that it can be written in terms 
of the four-dimensional anharmonic ratios [42] in an accordance with the results of refs. [38j. Thus it 
can be viewed as part of the remainder function, -R2->4, which is expected to depend only on the three 
anharmonic ratios Ui, i = 1, 2, 3. In section 3 we will come back for a closer look at this expression. 

As we have already mentioned in the introduction, this Mandelstam cut structure of six-point 
amplitudes in the multi-Regge kinematics can be generalized to a larger number of external legs, 
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Figure 5: The 2 — > 6 scattering amplitude. 



n > 6. This generalization exhibits the remarkable feature of integrability [15]. As an example, let us 
consider the eight-point amplitude in multi-Regge kinematics (Fig. [S]). Again, the amplitude 

can be written as a sum of terms which are compatible with the Steinmann relations. The number 
of terms is already 42 and will not be discussed here in further detail. We only mention that, for 
the fully-signatured amplitude, the Regge-pole contribution can also be rewritten in the factorized 
representation (cf . (|29"j) ) : 



A//" r lT2T3T4 



isipa v tit2 is 2 r 2 6 v t2ts is 3 r 3 



&^ T3T4 is 4 r£4 



(43) 



As it was already the case for the 2 — > 4 scattering amplitude, there exist Mandelstam cut contributions 
which appear only in special physical regions. The most interesting one is illustrated in Fig. [6l This 





Figure 6: The Regge pole and two- and three-gluon cuts in the 2 — > 6 amplitude. 



Regge-cut piece belongs to singularities in the angular momentum variables j'2, j'3, and j'4 of the t2, £3, 
and £4 channels, resp. Using the Mandelstam argument given above, it is easy to see that it appears, 
for example, in the physical region S3 > 0, S2345 > 0, s > 0, si < 0, S2 < 0, S4 < 0, S5 < [45]. This 
region is obtained by 'double twisting' of the planar amplitude and is further illustrated in Fig. [71 



1 for the 2 

^1=^2 = 1. 



n amplitude, the number of terms, 7V„, obeys the recursion relation N„ — YuZ=i NkN„-k with 
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Figure 7: The 'double twisting' of the 2 — > 6 amplitude. 

The detailed form of this contribution will be given in a forthcoming paper. Here we only mention 
that the singularities in j 2 and j'4 are described by the BFKL Hamiltonian of the two reggeon state in 
the adjoint representation, whereas the j'3 channel is governed by the spectrum of the BKP Hamiltonian 
of three reggeized gluons, projected on the adjoint representation. In section 3 we will give a more 
detailed discussion: in particular, it will be shown that this Hamiltonian is integrable and belongs to 
an open spin chain. 



3 BFKL approach and MHV amplitudes 

In this section we discuss six-particle amplitudes in the multi-Regge kinematics in some of the Man- 
delstam regions. We compare results obtained in the BFKL approach with those calculated using 
Wilson Loop/Scattering Amplitude duality at two loops. In particular, the analysis of the two-loop 
result allows to obtain the impact factor for the Mandelstam cut contribution beyond the LLA. We 
also discuss briefly the collinear limit and write an explicit analytic form of the all-loop remainder 
function in the Double Leading Logarithmic Approximation (DLLA). The section consists of two parts 
devoted to 2 — > 4 and 3—7-3 amplitudes. 



3.1 2-^4 amplitude 

The six-particle scattering amplitude corresponds to two physical processes, namely to 2 — > 4 and 
3 — > 3 scattering. Firstly we discuss the 2 — > 4 planar MHV amplitude illustrated in Fig. [8] and 
review the main result of the BFKL analysis applied to this case. The corresponding Mandelstam 
variables are defined as s = (pa +Pb) 2 , «i = (pa' + h) 2 , s 2 = (ki + k 2 ) 2 , s 3 = (p B > + k 2 ) 2 , S012 = 
(PA> + fa + k 2 ) 2 , S 123 = (p B > +ki + k 2 ) 2 , ti = (p A ~ PA') 2 , h = (p A ~ PA' ~ h) 2 , t 3 = (p B - PB') 2 
and the dual conformal cross ratios are given by 

ss 2 Slt 3 s 3 ti 
ui = , u 2 = — , u 3 = — . (44) 

■S012 S123 Soi2 12 S\ 23 t 2 

The multi-Regge kinematics, where 

- s > -S012, -s\ 23 > -si, -s 2 , -s 3 > -h, -t 2 , -t 3 > (45) 

implies 

l-m^+0, u 2 ^+0, n 3 ^+0, — ^— ~0(1), — ^— -O(l), (46) 

1 — U\ 1 — Ul 

which suggests that in this kinematics the convenient variables for the remainder function are 1 — u\ 
and the reduced cross ratios defined by 

u 2 „ u 3 

U2 = , u 3 = . (47) 

1 — Ul 1 — Ul 
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Figure 8: The 2 — > 4 gluon scattering amplitude. 



In the Regge limit they can be expressed through s 2 and the transverse momenta 

(k 1 + k 2 )2 x ^ k?q| a ^ klqf 



1 — U\ 



*2 



U 2 ~ 



(ki+k 2 )2q2 



2 ' 



u 3 ~ 



(ki+k 2 )2q2 



2 • 



(48) 



Note that « 2 and U3 are rational functions of cross ratios in four dimensions, but in the Regge limit 
they are simple cross ratios in the two-dimensional transverse space as one can see from (|48|) . The 
region of possible values of u 2 and U3 that correspond to physical momenta is depicted in Fig. |9] as a 
semi-infinite strip 



~ 2 

U 2 



~ 2 

U 3 



Figure 9: The physical values of \J~U2 and y[u% lie in the shaded semi-infinite strip. 

In the "Euclidean" kinematics, where all invariants are negative and thus all U{ are positive, the 
remainder function vanishes asymptotically as follows from the analysis presented in refs. j40| 142] . 
However, these studies also show that this is not the case in a slightly different physical region, where 
one or more dual confer mal cross ratios possess a phase. This happens when some energy invariants 
change the sign and here we consider one of such regions of the 2 — > 4 scattering amplitude having 

u\ = |ui|e _ * 27r , ii 2 and U3 are fixed and positive. (49) 

It corresponds to the physical region mentioned before (the Mandelstam region), where 

s,s 2 >0; si, s 3 , s i2) S123 < (50) 
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as illustrated in Figs. 0] and [TUJ The t- variables are all negative in the physical regions under consid- 
eration for the 2 — > 4 scattering amplitude in the Regge kinematics. It is worth emphasizing that the 
scattering amplitude in Fig. [Til] is still planar, but the produced particles have reversed momenta k\ 
and k 2 with a negative energy components. 




Figure 10: The Mandelstam channel of the 2 — > 4 gluon planar scattering amplitude. 

In the Mandelstam channel the remainder function grows with energy s 2 and was first calculated 
using the BFKL approach by two of the authors in collaboration with A. Sabio Vera in ref. |42j . 
The BFKL approach, based on the analyticity and unitarity was developed more than thirty years 
ago PP. In this approach one sums the contributions from the Feynman diagrams, which are enhanced 
by the logarithms of the energy (1 — u\ ~ (ki + k2) 2 /s2 in our case). The Leading Logarithmic 
Approximation (LLA) allows to write an integral representation of the remainder function R to any 
order of the parameter g 2 lns2- As it was already discussed in the previous section (see (|4"0|) ). the 
amplitude in this Mandelstam channel is given by [42] 

Af 2 ^ 4 = M^i R 2 ^ A = M*°l (1 + i A 2 ^ 4 ) , (51) 

where M^f is the BDS expression |35j and the correction A2->4 was calculated in all orders with a 
leading logarithmic accuracy using the solution to the BFKL equation in the adjoint representation. 
The all-order LLA expression for A 2 ^4 was given in (I4ip 



OO n QQ J 

- % E / -rV (™*r -t w +t - u ^ {u,n) - x ) • ( 52 ) 

Z n=-oo J -°° V + — 

Here fci, k 2 are complex transverse components of the gluon momenta, q\,q2,Q3 are the corresponding 
momenta of reggeons in the crossing channels. It is convenient to define holomorphic and antiholo- 
morphic variables in the transverse space as 



w = 1 j w = -p—- (53) 

k 2 qi k*q{ 

related to the reduced cross ratios of (f4T|) by 

I |2 ^2 U 2 \ i(<fo-6s) (A A \ !-^2-«3 l-Ui-n 2 -M 3 

\w \ = — = — , w = \w\e l(<P2 <P3; , cos(0 2 - 3 ) = ~ r^= - = n . • (54) 

us u s 2^u 2 u 3 2^u 2 u 3 

The energy behavior of the remainder function is determined by the Mandelstam cut intercept 

u(v,n) = -aE v>n , (55) 
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where a and E v n are the perturbation theory parameter and the eigenvalue of the BFKL Hamiltonian 
in the adjoint representation given by 

(X N [ 72 ( 72 \ / 72 \ 

° = ^T' K > n = ~2^^ +1>\l + i"+ VJ yj -2V(1). (56) 

Here ip(z) = T'(z)/T(z) and 7 = is the Euler constant. The two loop LLA expression for the 

remainder function in the BFKL approach was first found from (|5ip and (|52[) in ref. [42] and it reads 



R^^ LA = — ln(l — tii) ln£t2 In 113 = — ln(l — u\) In |1 + w\ 2 In 



i + I 

w 



(57) 



The remainder function in (|57p is pure imaginary and symmetric under w —> 1/w transformation, 
which corresponds to the target-projectile symmetry pa o pb, Pa 1 ^ Pb' and k± -f-)- fc 2 in accordance 
with (1521) . 

This result was shown by Schabinger [50] to agree numerically with the analytic continuation of 
the expression for the two-loop remainder function found by Drummond, Henn, Korchemsky and 
Sokatchev [51] from Wilson Loop/Scattering Amplitude duality. A rather complicated expression 
of ref. [5T] was largely simplified by Del Duca, Duhr and Smirnov [52 \ 153] and then by Goncharov, 
Spradlin, Vergu and Volovich (GSVV) [36]. The prediction in (|57p was analytically confirmed by two 
of the authors [48J performing the analytic continuation of the GSVV expression for the remainder 
function at two loops. The analytic continuation allowed also to extract the next-to-leading logarith- 
mic (NLLA) contribution, not yet available from the BFKL approach 

R mNLLA = in_ ln |w|2 ^2 (1 + w | 2 _ ™ ^3 1 1 + ^2 + ^ ln |^|2 ^_ w ) + L i 2 (-«;*)) 

-i27r(Li3(-tw)+Li3(-«;*)). (58) 

The NLLA remainder function in ()57[) is also pure imaginary and symmetric under w — > 1/w trans- 
formation. Both of the contributions are pure imaginary due to a cancellation of the real part coming 
from the Mandelstam cut, Regge pole and a phase present in the BDS amplitude as was shown by one 
of the authors |47] , Starting at three loops the cancellation does not happen anymore and the real 
part gives a non-vanishing contribution at the next-to-leading level. The analysis of ref. |47] based 
on analyticity and other general properties of the scattering amplitudes resulted in a formulation of 
the dispersion-like relation for the real and imaginary parts of the remainder function in the Regge 
kinematics in this Mandelstam region 

/ioo j 
-/Me-™(l-«ir, (59) 
-ioo 

where the first term in RHS corresponds to the contribution of the Regge pole (see (|34|) ). This term 
as well as the phase 5 in LHS of (|59p are obtained directly from the BDS formula (see also (|36|) ) 

5 = — Hu 2 u 3 ) = —ln ]TT - ¥ , Wa6 = — ln^ = T lnH . (60) 

The second terms in RHS of (|59p stands for the contribution of the Mandelstam cut (see (|35fl ). The 
coefficient 7^ ~ 4a is the cusp anomalous dimension known to an arbitrary order of the perturbation 
theory. The only unknown piece in (|59p is the real function f(uj), which contains the Mandelstam cut 
in oj, depends only on the transverse particle momenta and has no energy dependence. In the leading 
logarithmic approximation f(co) is given by 

f LLA [u) * y r du 1 (^)--l M -+f (6 i) 
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where u{y,n) is defined in (|55p . 

The dispersion-like relation in (|59p was used [Hj] for calculating the three loop contributions to 
-^2->4 (leading imaginary and the sub-leading real terms) in the multi-Regge kinematics 

E!i\ L A LA = iA^Ja 3 = ^ ln 2 (l - « a ) (in \w\ 2 In 2 |1 + - | In 3 |1 + w\ 2 (62) 
-^ln 2 | w | 2 In 1 1 + w\ 2 + i In |w| 2 (Li 2 (-w) + Li 2 (-«;*)) - Li 3 (-u;) - Li 3 (-u;*)^ 

and 

K (4 3 -if LLA ) = y ln ( X - n i) ( ln l^l 2 ln ' I 1 + - I In 3 |1 + w\ 2 (63) 

~ In 2 \w\ 2 ln |1 + w\ 2 - \n\w\ 2 (Li 2 (-w) + U 2 (-w*)) + 2Li 3 (-w) + 2Li 3 (-u>*)^ . 

As in the two loop case, both (|62p and (|63p are symmetric under w —> 1/w transformation, which 
is obvious from (|52p and corresponds to the target-projectile symmetry of the scattering amplitude. 
The corrections, subleading in the logarithm of the energy, are not captured by (|52|) and require some 
knowledge of the next-to- leading impact factor and the eigenvalue of the BFKL Kernel in the adjoint 
representation. 

While the latter is still to be found from the next-to-leading BFKL equation constructed by Fadin 
and Fiore [541 155] . the correction to the impact factor was obtained in ref. |49] extracting it from (j58|> . 
The integrals in (|52p and (|59p come as a convolution of the propagator of the BFKL state Gbfkl 




Figure 11: A graphic representation of the expression in (I52p . Two impact factors xi and X2 are 
convoluted with the propagator of the BFKL state Gbfkl- 

and two impact factors xi and X2 as shown in Fig. [TTJ The leading logarithmic impact factor xf LA 
was calculated by two of the authors directly from the Feynman diagrams in ref. 



x\ LA 



1 



1 



2 (iu + 5 



si 



-iv+i 



LLA 
< A2 



1 



1 



2(w/-a) \k* 2 



while the NLO impact factor was extracted from (i58j) and read [19 

/ 



xf LO 



E, 



\ 



n 



LLA 
Xl , 



(64) 



(65) 
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where E v n is defined in (|56|) . The NLO correction to xi has a similar form found in ref. [49 j. An 
important feature of x^ LO in f|65[) is the fact that, in contrast to the leading order, it has lost the 
property of holomorphic separability: we cannot write the factor in front of X\ LA i n dS5J) as a sum 
of terms, which depends only on either iv + n/2 or —iv + re/2. It is worth emphasizing that the 
NLO impact factors xf are factorized in the product of the Born impact factors in (|64[) and a term 
expressed through the eigenvalue E v n of the BFKL equation in the LLA. The form of the NLO impact 
factor in the u, n representation resembles the three- loop remainder function in the LLA, emphasizing 
the intimate relation between the two. Indeed, it is easy to see that expanding the integrand of (|52p 
to the third order in a one gets the E 2 n term. 

In the general case the integral in (|52p is not easy to calculate, but one can consider a more 
restrictive kinematics, where it can be found explicitly at any order of the coupling a. One of such 
possibilities is the so-called collinear kinematics, when two adjacent particles become collinear, e.g. if 
in Fig. [8] the momenta ps and Pb 1 coincide. In the limit t% — > the remainder function vanishes at two 
loops and beyond, in both the direct channel of Fig. [8] and in the Mandelstam channel of Fig. [TU1 The 
multi-Regge limit followed by the collinear limit in terms of the dual conformal cross ratios (compare 
to the Regge kinematics in (j46j) ) reads 

1 - ui -»■ +0, u 2 -»■ +0, u 3 -»■ +0, — ^— = u 2 -»■ +0, — — = « 3 ~ 1, (66) 

1 — U\ 1 — U\ 

which in terms of w and w* implies (see (|54j) ) 

l-m^+0, H-»-+0, cos(<p 2 -0 3 ) - 0(1). (67) 

For \w\ — > the main contribution to the LLA remainder function given by (|52p comes from poles 
at v = ±m/2 for the conformal spin re = 1. In this case one can drop the ip functions in (|56p and 
perform integration of (|52p at any order of the coupling constant in the Double Leading Logarithmic 
Approximation (DLLA), where one sums contributions from the powers of a In \w\ ln(l — u\) resulting 
in [56] 

R^ A = 1 + i27racos(^ 2 - </> 3 ) \w\ (l - I foy/ala \w\ ln(l - uM , (68) 

where I n (z) is the modified Bessel function. The Double Leading Logarithmic Approximation is 
analogous to the summation of the contributions of powers of g 2 In s In Q 2 in Deep Inelastic Scatter- 
ing (DIS), where the BFKL and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations 
overlap. Then the real part of the next-to-leading corrections to the remainder function in DLLA was 
calculated [39] using the dispersion relation ([59]) and ([6"8"|h 



3.2 3—7-3 scattering amplitude 

Beyond the 2 — > 4 scattering the six-particle scattering amplitude describes also the 3 — > 3 scattering 
illustrated in Fig. [12j The analysis of the 3 — > 3 amplitude is similar to that of the 2 — > 4, but there 
are some rather interesting features we want to emphasize. Firstly we start with the definition of the 
kinematic invariants si 3 = (ps + &i) 2 , S02 = (lM' + &2) 2 ,s = (pb + h +pa) 2 , t 2 = (pa ~PA' ~ k 2 ) 2 , si = 
(A;i +p A ) 2 ,s 3 = (p B > + k 2 ) 2 ,t 2 = (p A - Pa' + h) 2 ,h = (PA - PA') 2 and t 3 = (p B - Pb 1 ) 2 ■ The dual 
conformal cross ratios are expressed in terms of these invariants as follows 

•S13S02 t\H S\S 3 

ui = — -;— , u 2 = —r, u 3 = — — . (69) 
s t 2 t 2 t 2 s t 2 

In the multi-Regge kinematics for the direct channel in Fig. [T2l where all invariants are negative 

- s » -a x , -s 3 , -t 2 > -tx, -t 2 , -t 3 > 0, (70) 
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Figure 12: The 3 — > 3 gluon scattering amplitude. 



the remainder function R^^ is zero, while in the physical region of the Mandelstam channel depicted 
in Fig. [T3l where 

si,s 3 ,si3,s 2 < and s,t' 2 >0, (71) 
it contains a non- vanishing contribution, growing with energy t 2 - In the Mandelstam channel (I7ip in 




Figure 13: The 3—^3 gluon scattering amplitude in the Mandelstam channel given by s\, S3, S13, S02 < 
and s, t 2 > 0. 



the multi-Regge kinematics the dual conformal cross ratios (169j) possess a non-zero phase 

ui -)• \ui\e t2n , u 2 ->■ |w2|e 17r , U3 ->■ |xi 3 |e i7r 
and the analytic continuation of the GSVV expression in multi-Regge kinematics gives 



R 



(2) LLA+NLLA 
3^3 



ITT 



\n(ui — 1) In 1 1 + w\ 2 In 



1 + 



1 



w 



2 vr 2 



+ — In II + w\ 2 ln 
2 1 1 



1 + 



w 



(72) 



(73) 



t7T ITT 

-— ln|w| 2 ln 2 |1 + w\ 2 + — In 3 |1 + w| 2 - iirln\w\ 2 (Li 2 (-w) + Li 2 (-u>*)) 
+i2ir (Lia(-u;) + Li 3 (-w*)) . 

As in the 2—7-4 case the remainder function has the target-projectile symmetry (pa <-> Pb, Pa' ^ Pb', 
k\ ■<->■ ^2 or I «; I — » l/|to|), but in contrast to the 2 — )• 4 amplitude ([73]) has a real part In 1 1 + 
In 1 1 + -M 2 . This fact is in full agreement with the dispersion-like relation |47| for the 3 — > 3 



tr 



amplitude 



R 



3^3e 



COS lTU) a b 



doj 

2^1 



(74) 
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which differs from (|59p by signs of the phase on LHS and the integral on RHS, as well the absence of 
the phase e~ l7TUJ in the integrand. This phase mixes the contribution from the Regge poles and the 
Mandelstam cut in the dispersion relation (|59|) , which leads to the full cancellation of the real part at 
two loops for the 2 — > 4 amplitude. In the 3 — > 3 case this cancellation does not happen anymore due 
the absence of the phase e~ mu , and one obtains a real term in (|73p . The real part does not cancel out 
in both the 2 — > 4 and 3 — > 3 amplitudes at higher loops. The dispersion relation (|74p was used by 
two of the authors [49] to find the LLA and the real part of the NLLA contributions to the remainder 
function of the 3 — > 3 amplitude 

R ^LLA = _™ ln2( ^ _ 1} / H2 ^2 |x + _ 2 ^3 ^ + ^ ^ 
1 1 

+-ln|u>| 2 (Li 2 (-u>) + Li 2 (-w*)) - - In 2 \w\ 2 In |1 + w\ 2 -Li 3 (-iy) -Li 3 (-tt>*) 

and 



* (4 3 _) 3 ™) = In( ttl - 1) (in 2 |1 + w\ 2 In 



1 

1 + - 

w 



+ ln|l + u;| 2 ln 2 



1 




1 + - 




w 





(76) 



As we have already mentioned in the 3 — > 3 case there is no mixing between the real contributions 
coming from the Regge poles and Mandelstam cuts. This fact allows us to make a prediction valid at 
an arbitrary value of the coupling constant for the following expression [57] 



Vt[R 3 ^ 3 e- l7TS ) =cosiroj ab . (77) 



Both S and uj a b are known and are given by (|60p as functions of the cusp anomalous dimension and the 
anharmonic ratios. However there is a difficulty in understanding (I77p at the strong coupling because 
of the rapid oscillations as a — > 00. 

Similarly to the 2 — > 4 case one can also consider the collinear limit t% — > (jw\ — > 0) preceded by 
the Regge limit. In this kinematics it is possible to calculate explicitly the LLA and the real part of 
the NLLA 3—7-3 remainder function at an arbitrary number of loops [56J . 

In the next section we return to the discussion of section [2] and consider a composite state of an 
arbitrary number of reggeized gluons. These BKP states appear in the scattering amplitudes with 8 
or more external legs. 



4 Integrability of the n gluon Hamiltonian 

Here we discuss composite states of n reggeized gluons in the adjoint representation at large N c (cf. a 
similar approach for the simple case n = 2 in ref. |42j). One can write the homogeneous BKP equation 
for its wave function described by an amplitude with amputated propagators in the form [45] 



H* = EV, A n = -^E. (78) 

Here H is a redefined hamiltonian obtained after a subtraction of the gluon Regge trajectory cj(t) 
containing infrared divergencies. Namely, the Regge trajectory of the composite state is |45] 

u n {t) = a Q - In +A n , a = , (79) 

where A n is the infrared stable quantity expressed in terms of the energy E. 
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The hamiltonian H in the multi-color limit can be written in the holomorphically separable form 
(see 05]) 

n— 1 n 

H = h + h* , h = In 1 2 W + ^ /i*,r+i > 9 = ^ , (80) 

^ r=l 1 

where the pair hamiltonian r+l is transposed to the corresponding unamputated operator (fT7|) 

/i* r+1 = ln(p r p r+ i) + p r ln(p r r+ i) — + p r+ i ln(p r r+1 ) — h 27 . (81) 

Pr Pr+1 



It is seen from (|80p that the holomorphic hamiltonian for the composite state in the adjoint represen- 
tation differs from the corresponding expression for the singlet case h\ ' (|17p after its transposition 
only by the substitution 

h n>1 ^]n^, (82) 
Q 

which is related to the fact, that the planar Feynman diagrams have the topology of a strip and the 
infrared divergencies in the Regge trajectories of the particles 1 and n are not compensated by the 
contribution from the pair potential energy V n \. 

It turns out, that the eigenvalues E do not depend on |g| 2 due to the scale invariance of H, as it 
will be demonstrated below. As a result, the t-dependence of u; n (t) is the same as in the gluon Regge 
trajectory. 

The normalization condition for the wave function in two-dimensional space can be written as 
follows 

/n—l n n 

Hd 2 Pr y*ll\p s \- 2 y, J2 Ps = q. (83) 
r=l s=l s=l 

Using the duality transformation (see [9] and [45] ) 

d 

Pi = ^0,1 , Pr = Z r -l,r , q = ZQ,n , Pr,r+1 = « = id r , (84) 

the holomorphic hamiltonian can be presented as follows 

n-i 

h = ^ Z -^^ + J2K >r+1 . (85) 

z 0,n r= i 

Further, by regrouping its terms we can write the holomorphic hamiltonian in another form [45] 

n-2 

h = -21nz ,n + H4,l 9 l) + ln (4-l,rA-l) + 27 + Kr+l > ( 86 ) 

r=l 

where the new pair hamiltonian is 

K tr+ i = ln(z 2 r+1 d r ) + ln(^ )r+1 (9 r+ i) - 21nz r)r+ i + 27 

= ln(d r ) + ln(«9 r+ i) + — - In z r , r+ i d r + In z r . ir+1 d r+ i + 27 . (87) 

Or Or+1 

The operator h' r r+1 coincides in fact after the substitution z r — > p r with the corresponding hamiltonian 
in the coordinate representation (I17h acting on the wave function with non-amputated propagators. 
It is important, that h (j86[) is invariant under the Mdbius transformations 

*k - (88) 
czfc + a 
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and does not contain the derivatives 8q and d n . Therefore we can put 

z = , z n = oo , (89) 
which leads to the simplified expression for h 



n-2 



h^ti = ln(zf5i) + In(SW_i) + 27 + ^, r +i • (90) 

r=l 

To return to initial variables in the final expression for the wave function one should perform the 
following substitution of Zk 

Zk~ z n q - ^ r=1 Pr 

According to the above representation ([831) f° r its transposed part h' 1 can be obtained from h 
by the similarity transformation which can be written in terms of b! as follows 



'n-2 \ _1 /n-2 \ 

i n^+M ^ zi ( n^+i) 1 ( 92 ) 

which is compatible with the following normalization condition for the wave function in the full two- 
dimensional space 

ii*ii?=/wnii%F |<f|2 - (93) 

J |~1 1 ) .y \^r,r+l\ 



On the other hand, from the expression (|90p for /i' we obtain another relation for h 

/n—l \ /n—l \ ~ 



n m m n M ' (94) 



vr=l / \r=l 



corresponding to the second normalization condition for ^ compatible with the hermicity properties 
of the total hamiltonian 



, n—l n—l 



11*111 = /"n^** n i^i 2 *- 

J r=l r=l 

By comparing two above relations between b! and h n one can conclude (cf. |1U|). that the operator 

n— 2 n—l 

A' = *1 J] Z SjS+ i J] d r (96) 

s=l r=l 

commutes with the holomorphic hamiltonian 

[A',h']=0. (97) 

4.1 Integrable open spin chain 

Let us verify |45| . that the holomorphic hamiltonian b! (|90p also commutes with the differential 
operator D(u) being the matrix element T22 of the monodromy matrix (cf. |10j ) 

T(u) = ( ^ ) = L x (u)L a (u)...Z te -i(«) , (98) 
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where the L-operator is defined by the relation 



Lr(u)=[ U +Sf r J d L* ) • <^> 



-iz^d r u — iz r d r 

To prove the commutativity of h! and D(u) one can use the following relation 

[L k (u) L k+ i(u), h' kjk+1 ) = -i {L k (u) - L k+ x(u)) , (100) 

valid in particular due to the Mobius symmetry of the pair hamiltonian 

[M k , k +i, h' k>k+l ] = , M Kk+1 = M k + M k+1 . (101) 

and the fact that its eigenvalue is a linear combination of polygamma functions (see (|13p ) 
Relation (|100|) leads to the equality 

n-2 

[T(u),y^ y h' rj . +1 ] = iL 2 (u)L 3 (u)...L n -i(u) - iLi(u)L 2 (u)...L n _ 2 (u) • (102) 

r=l 

On the other hand, one can easily verify, that 

[T 22 (u), ln(4di)+ In d n ^] = (0, 1) [T(u), Hzfdi) + In d n _i] ( ° 

= -i(0, l)(L 2 (u)L 3 (u)..X n - 1 (u)-L 1 (u)L 2 (u)...L n - 2 (u)) ( ? ) , (10:5) 



x 1 

which proves that the differential operator D(u) = T 22 (u) is an integral of motion 

[D(u),ti] =0. (104) 
Note, that, if instead of condition (|89p we shall use the equivalent condition 

z = oo , z n = , 

the matrix element ^4(u) = Tn(u) of the monodromy matrix will be an integral of motion. 

Thus, our hamiltonian is the local hamiltonian for an open integrable Heisenberg spin model with 
the spins which are generators of the Mobius group. 

With the use of the following decomposition of the L-operators 

Lr{u)=( I ) + ( 1 ){Zr, 1) id r (105) 



u J \ —z, 

one can construct the matrix element T 22 = D(u) in an explicit way 

n-1 

D(u) = Y J U n - 1 - k q' k , (106) 
k=0 

where 

fc-l k 

Qk = ~ z nYl z rs,r s+ i Ylidn ■ (107) 

0<r±<r2<...<r k <n s=l t=l 

Note, that one can parameterize the monodromy matrix in another form 

= { Uu) Jo{u) - Mu) ) - = ± • ( 108 ) 

In this case the Yang-Baxter equations for the currents jn have the Lorentz-invariant representation [9] 

lUu),j u (v)} = {f(u)f{v) ~ f{v)f{u)) . (109) 

Here e^pa is the antisymmetric tensor in the four-dimensional Minkowski space and €1230 = 1 , 9p,v = 
(1,-1,-1,-1). 
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4.2 Composite states of two and three gluons 

In the case n = 2 we have only one non-trivial integral of motion 

4 x = -%z x d x . (110) 
Taking into account the normalization condition for the eigenfunction in the two-dimensional space 

= |^|t|f, (111) 

we find the orthonormalized and complete set of eigenfunctions 

T (2) -i+m / 1+n. _ 1+n. 

& m L= Zl ' {zt) * +m , m=—^-+tu, m =——-iv, (112) 

satisfying the single- valuedness requirement. Note, that using the substitution (f91~j) one can reproduce 
the wave functions of two gluon composite states in the momentum space (see [42J). 
For the case n = 3 the operator D(u) is given below 

D 3 (u) = u 2 - iu(zi di + z 2 d 2 ) + Z\.z\$. d\d 2 ■ (113) 

With taking into account the normalization condition 



■J- 1 1 2 f d 2 zi d 2 z 2 . _ l2 



|Zir|2l,2'° 

one can search the holomorphic eigenfunction of this operator in the form 



i 



^ = V +m /(f ) • (115) 



The function f(x) satisfies the equation 

x(l - x)d 2 + (- + m)(l - x)d + A ) / = , x = — , (116) 
2 / zi 

where A is the eigenvalue of the operator z\Zi^ 2 d\d 2 . Two independent solutions of this equation can 
be expressed in terms of the hyper geometric function F 

flips) = F( ai ,a 2 ; 1 + ai + a 2 ; x) , f 2 (x) = x ai+a2 F(-a 2 , -au 1 - a x - a 2 ; x) , (117) 

where the parameters a\ and a 2 are obtained from the set of equations 

1 

«i + a 2 = — - + m , a\a 2 = -A . (118) 

We have the similar solutions for the eigenvalue of the operator D* in the antiholomorphic subspace. 
They can be obtained by the substitution 

1 — n 

x — > x* , a± — >■ Hi , a 2 — > a 2 , m — > m = — h iv . (119) 

To construct the wave function * with the property of the single-valuedness in the two-dimensional 
subspaces z\ and x we should write a bilinear combination of the functions fi(x) and the corresponding 
functions in the anti- holomorphic subspace fi{x*). 

One can write the integral representation for the wave function satisfying the above constraints 

* ~ (z*r^ J ^ y --(y*)-- (^)" (f^Y ,* = I , (120) 

where the integration is performed over the two-dimensional plane y. Note, that the integrand has no 
ambiguity in the points y = 0, l,x due to the additional constraints for the parameters a^cTj [45J 

ai - a! = iV ai , a 2 -a 2 = N a2 , (121) 

where A r ai , A^ a2 are integers. Moreover, the function * near the points x = 0, l,oo can be presented 
in terms of the sum of products of hypergeometric functions in (|117p . 
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4.3 Hamiltonian and integrals of motion 

The holomorphic hamiltonian for composite states of two reggeized gluons can be written as follows 

h = Hzldx) + ln(di) + 2 7 = %l>(z x dx) + i>{-z x dx) + 2 7 . (122) 

Acting by h on the function zf we obtain 

hz{ = e{5) 4 , e(5) = if,(S) + tf>(-8) + 2 7 . (123) 

In the case of wave function ()112p satisfying the single-valuedness and orthonormality conditions in 
the two-dimensional space one derives the following expression for the total energy |42| 

E m ,m = e m + £m, e m = + m) + - m) + 2 7 . (124) 

Note, that it does not coincide with the corresponding result (I13p for the Pomeron state. 

For the case of composite states of n reggeized gluons the holomorphic hamiltonian (|86p in the 
region 

zi < z 2 < z 3 < ... <C . (125) 
becomes the sum of the disconnected hamiltonians 

n-1 

h' = J2 (^Mi) + i>{-z r d x ) + 2 7 ) . (126) 

r=l 

As a result, we obtain, that the wave function in this limit is factorized |45j 

n-1 

*a 1 ,a 2 ,...,a„_ 1 = n<" ( 127 ) 
r=l 

and the energy is the sum of the particle energies 

n-1 



e = 

r=l 



^e(a r ). (128) 



The parameters a r for these solutions and a r for anti-holomorphic solutions are obtained from the 
single-valuedness condition and the normalizability 



where f r are real and n r are integer numbers. 

The eigenvalue of the integral of motion D(u) is expressed also in terms of these parameters 



n-1 



A(«) = Jj(tt-ior). (129) 

r=l 

4.4 The Baxter-Sklyanin approach 

To find a solution of the Yang-Baxter equation for the open spin chain one can use the Bethe ansatz. 
For this purpose it is convenient to work in the transposed representation for the monodromy matrix 
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where the L-operator can be chosen as follows 



1 1 , , I u -\- id r z r id r i , 

L r(«) = .-a J2 .. jo „ • (131) 



-id r Zr u — id r z. 
The pseudo- vacuum state is defined as a solution of the equation 

= 0. (132) 

It can be written in the form [12] 

n-1 

= J] ^ • ( 133 ) 

r=l 

Note, that the function l^ol 2 does not belong to the principal series of the unitary representations. 
As a result, the states constructed in the framework of the Bethe ansatz by applying the product of 
the operators j r + {u r ) to 

k 

*fc = IIi+(«r)*o (134) 

r=l 

are non-physical. Nevertheless, these states are eigenfunctions of the integral of motion 

D\u)¥ k = (j*(u) - Mu))*i = A(«)*| (135) 

providing that 

k 

A(„) = (. + ,)-' II _L- S (, + (136) 

is a polynomial, which leads to a quantization condition for the Bethe roots ut- If we parameterize 
this polynomial as follows 

n-l 

A(«) = JJCti-io,), (137) 

2=1 

the above defined Baxter function Q(u) can be calculated 



ax r(-t« + i) 

As it was mentioned above, the polynomial solutions for Q(u) are non-physical, because the cor- 
responding wave functions \& do not belong to the principal series of unitary representations of the 
Mobius group. We should find a set of non-polynomial solutions Q s {u) satisfying this physical re- 
quirement. 

According to E. Sklyanin |15j the correct variables in which the dynamics of the Heisenberg spin 
model is drastically simplified are the zeroes b r of the operator B(u) = j+(u) entering in the mon- 
odromy matrix 

n— 2 ra— 1 

B{u) = P n _! \\ (U - b r ) , P n _! = % r , (139) 
k=l r=l 

where the operators b r and P n -± commute each with others 

[6 P ,6 a ] = [6 r ,i' n _ 1 ] = 0. (140) 

It is convenient to pass from the coordinate representation z to the Baxter-Sklyanin representa- 
tion [13J , in which the currents j\(u) and {j\(u))* (together with the operators b r , b* and P n —i, Pn-i) 
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are diagonal. We denote the eigenvalues of the Sklyanin operators by b r ,b*. The kernel of the uni- 
tary transformation to the Baxter-Sklyanin representation is known explicitly for the cases n = 2, 
n = 3 and n = 4 [13] . For general n this integral operator can be presented as a multi-dimensional 
integral [14J. 

In the Baxter-Sklyanin representation the wave function in the holomorphic subspace can be 
expressed as a product of the pseudo-vacuum state in this representation ^>o(-fn-i ; bi, b 2 , — > b n - 2 ) and 
the Baxter functions Q(ut) 



n-2 

n — 1 



**(P n _ 1 ;6 1 ,...,6 n _ 2 ) = P n _ 1 2 m nQ(6 fe )M/ (P„-i,&i,...A- 2 ), (141) 



fc=i 



where the power- like behavior in the variable P n -\ is in an agreement with the normalization condition. 
The analogous representation is valid for the total wave function 



n-2 



vP + (P n _ i; h, b n _ 2 ) = P n _l m (P*^)- 2 ^- m [J Q&) *o(Pn-i; &i, 6n- 2 ) (142) 

fc=i 

with the use of the generalized Baxter function Q(u) being a bilinear combination of the usual Baxter 
functions in the holomorphic and anti-holomorphic subspaces 

Q(u)=Y / d s ,tQs(u)Qt(u*). (143) 

s,t 

Here Q s (u) are different solutions of the Baxter equation with the same eigenvalue A(u). The coeffi- 
cients d Si t are chosen from the requirement, that the function Q(u) satisfies the normalization condition 
everywhere including the points where the functions Q s (u) and Qt(u*) have the poles |13|. 114] . For 
the periodic spin chain this condition leads to the quantization of the eigenvalue of the operator 
A(u) + B(u) although a simpler method of quantization is based on the requirement, that all Baxter 
functions corresponding to the same eigenvalue should have the same holomorphic energy [13]. In the 
case of the open Heisenberg spin model the situation is simpler and will be discussed below. 

4.5 Baxter-Sklyanin representation for three gluon states 

For the states composed from three reggeized gluons the transposed integral of motion in the holo- 
morphic subspace is 

D\{u) = u 2 - iuidx z x + d 2 z 2 ) + d x d 2 z x z x ,2 (144) 
and the operator f\ is given below 

f + = iu(di + d 2 ) - d x d 2 z l2 = %{dx + d 2 ) (u - Si) , (145) 

where 

di + d 2 

The operator j l + is easily diagonalized after a transition to the momentum representation, where 

idl fpi,p 2 = Pi fpi,P2 > ^2 fpi,p 2 = P2 fpi,p 2 ■ (147) 

In this case the eigenvalue equation for jt_ has the form 

d d \ 

u(pi +P2) - ipW2{t, -k— ) / = (pi +P2)(u - bi) f , (148) 

dpi op 2 J 
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where b\ is the eigenvalue of b\. Its solution is given below 

f = x(P!+P2,h) 0£\ \ (149) 

where \ is an arbitrary function of p\ + p 2 and bi. The dependence of from p\ +P2 is fixed by the 
normalization condition 

(Pi + P2)~ ai - a2 ■ (150) 

On the other hand, the eigenvalue equation for the integral of motion in the momentum space can 
be written in the form 

d ( d d \ 

PlP2 dp~i \dp~2 ~dp~J * (P1,P2) = aia2 ^ (P1,P2) • (151) 

Using the ansatz 

*(pi,p 2 ) = (pi + P2)- ai ' a2 v(y) , y = -, (152) 

pi 

we obtain the following equation for the function rj(y) 

(y 2 3 2 + { ai + a 2 + l)yd- ai a 2 ) v(y) = {~V 3 d 2 -2y 2 8) V (y) . (153) 
There are two independent solutions of this equation [15] 

)) i(|/) = -F(l-ai ) l-a 2 ,2;--) (154) 

y y 

and 

, . r(-ai)r(+a 2 ) „ „, i , . 

m(y) = - rn , [ y~ ai F(- ai , l- ai ,l + a 2 - ai ;-y). 155 

T(l + a 2 - ai) 

One can construct the bilinear combination of these solutions having the single-valuedness property 
in the y-space. Finally with the use of the integral representation for the hyper geometric function the 
wave function ^ in the momentum space can be written as follows 



*\puP2) = (pi +P2)- ai - a2 {p\+Pir~ ai -~ a2 4>{y) , (156) 

where 4>(y) is given below 

m = J d 2 t (^+i) ai (^ + 1 ) ai t 1 -^ 2 " 1 ^-**)" 2 " 1 ( 157 ) 

and satisfies the single valuedness condition in the y-space due to the quantization conditions for a r 
and a r . 

The transition to the Baxter-Sklyanin representation (u, u) corresponds to the Mellin-type trans- 
formation of 4>(y) 

cf>(u, u) = |0 y~ iu (y*)" l >(y) . (158) 
The inverse transformation corresponds to the Baxter-Sklyanin representation for the wave function 

—iu / * \ —iu 



where 



// \ — IU / * \ — iu 

d 2 u<f>(u,u) (JiJ i^J , (159) 

N N f f°° 
-iu = iv u + ~y , -iu = iv u ^ , J d 2 u = J dv u . (160) 



N u =— oo 
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One can interpret the wave function 4>(u, n) in the Baxter-Sklyanin representation as a product of 
the pseudo- vacuum state nn and the total Baxter function [35] 



4>(u,u) = uuQ(u,u) , (161) 

where 

T(iu)T(iu) T(-iu- a^Tj-iu- a 2 ) . , 

WK ' ' ~ r(i - »«) r(i - iu) r(i + m + s x )r(i + «z + s 2 ) ' 1 J 

This expression can be written in the factorized form 

Q(u,u) ~ Q(u,ai,a 2 ) Q(u,ai,a 2 ) , (163) 

where 

^, . r(— in — ai)r(— iu — ao) . < , 

0(«,ai,a 2 =— r2 /l ■ n 164 

1 ^(1 — zn) 



$( u ) = J»"H'H-»"- «y; -v;r^ _ ^ (165) 



/sin(7r(— in — ai)) sin(7r(— in — a 2 )) 
sin 2 (— inn) 

The expression Q(u, a\, a 2 ) differs from the Baxter function in the holomorphic space 



n( \ r(-w-oi)r(-w-a 2 ) , , 

= f^T^ (166) 

only by the periodic function <3?(n) and therefore it can be considered also as a Baxter function. This 
additional factor can be included in the definition of a new pseudo- vacuum state. Really this pseudo- 
vacuum state can be considered as the additional factor for the wave function in the Baxter-Sklyanin 
representation providing correct hermicity properties of the hamiltonian and integrals of motion in 
this representation!! (see also ref. |14|). 



5 Conclusion 

In this review article we have outlined the role of Mandelstam-cut contributions in the remainder 
functions for the BDS amplitudes. Particular emphasis has been given to the integrability of the 
Hamiltonian which describes the energy spectrum of the states of n reggeized gluons. These cut 
contributions appear in multi-Regge kinematics in special physical regions, where some energies are 
negative. For the cut corresponding to the composite states of n reggeized gluons the number of ex- 
ternal particles should be k > 2 + 2n. The wave functions of these states in the adjoint representation 
satisfy the BFKL-like equations, which have the property of holomorphic factorization and are inte- 
grable in LLA. The corresponding holomorphic hamiltonian coincides with the local hamiltonian for 
an integrable open Heisenberg spin model. The Baxter equation for this model is reduced to a simple 
recurrence relation and can be solved in terms of the product of the T-functions. We constructed the 
wave functions of composite states of 2 and 3 gluons explicitly. 

We thank L. D. Faddeev, V. S. Fadin, E. M. Levin, J. Maldacena, A. Sabio Vera, V. Schomerus, 
A. Sever, M. Spradlin, Chung-I Tan, C. Vergu, P. Vieira and A. Volovich for helpful discussions. 

2 We thank Prof. F. Smirnov for discussions related to this important interpretation of the pseudo-vacuum state. 
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